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Abstract. We present some results obtained from the synthesis of Stokes 
profiles in small-scale flux tubes with propagating MHD waves. To that aim, re- 
alistic flux tubes showing internal structure have been excited with 5 min period 
drivers, allowing non-linear waves to propagate inside the magnetic structure. 
The observational signatures of these waves in Stokes profiles of several spectral 
lines that are commonly used in spectropolarimetric measurements are discussed. 



1. Introduction 



The observed wave behaviour in network and facular solar regions, and the role 
of waves to connect photospheric and chromospheric layers, have drawn a con- 
siderable discussion in the recent literature, since it has been proposed that thin 
magnetic flux tubes in these regions can act as wave guides, supplying energy 
to the upper layers, in particular for waves with periodicity in the 5-min range 



( De Pontieu et al1l2004l ). However, observational studies of the relation between 



the photospheric and chromospheric signals do not provide a unique picture. It 
i s not clear whet h er the oscillations r emain coherent through the atmospher e 
fiLites et all Il993l : iKriiger et all l200ll : lJudge et all l200ll : ICenteno et aP l2006l ). 
There are some hints that the couplin g between the pho t osphere and the chro - 



mosphere is within vertical channels (Judge et alJ l2Q0ll : ICenteno et al. 2006lV 
howeve r, other results point toward inclined wave propagation (IBloomfield et al 
I2OO6I : iDe Pontieu et al.ll2003l ). In particular, all these studies point out the im- 
portant role of the magnetic canopy and the cs = va transformation layer. The 
dominance of long-period oscillations in the chromosphere of these regions, also 
remains to be explained. 

Here we perform numerical simulations of magneto-acoustic waves in small- 
scale flux tubes, with a subsequent Stokes diagnostics. We address the ques- 
tions of coherence of oscillations through the atmosphere, wave behaviour at the 
canopies and the change of the wave period with height. 



2. Simulations and spectral synthesis 



We have solved n umerically the set of non-linear 2D MHD equations using the 
code described in iKhomenko fc ColladosI (|2006l . [20071 ). As initial condition, we 
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Figure 1. Snapshots of the vertical velocity (normalized to ) in adiabatic 
simulations excited with a vertical driver (left) and a horizontal driver (right) 
after about 800 sec elapsed time. The vertical lines represent the magnetic 
field lines. The two horizontal lines are contours of constant c^/v\^ the thick 
line corresponding to va = cs and the thin line to c|/v^ = 0.1. 



have used a magnetost atic flux tube model constructed following the method by 
Pneuman et all (| 19861 ). with a radius of 100 km at the base of the photosphere. 
This model is a second order approximation to a thin flux tube and has both 
vertical and horizontal variations in all atmospheric parameters. The magnetic 
field along the flux tube axis varies with height from 740 G in the photosphere 
to 37 G at z=2000 km, and the plasma (3 decreases from 4 to 7xl0~^. We 
compare here three simulation runs. In the first one, oscillations were excited by 
a photospheric driver with horizontal velocity given by Vx{x^ t) — Vq sin(c<;t) x 
exp(-(x - xo)^/2a^), where a = 160 km, 27v/uj = 300 sec, Vq = 200 m s-^ 
is the initial velocity amplitude at the photosphere and xq corresponds to the 
position of the tube axis. This run was performed in an adiabatic regime. In 
the other two runs, oscillations were excited by a periodic photospheric driver in 
vertical velocity with the same horizontal and temporal dependence as before. 
One of these runs was adiabatic while the other was carried out with a radiative 
relaxation time of trr = 10 s (constant through the whole atmosphere). The 
radiative losses were described by Newton's law of cooling. The Stokes spectra of 
the Fe l pair of lines at 1.565 fim a nd of the Si I line at 1.083 /im were sy nthesized 
with the help of the SIR code fjRuiz Cobo fc del Toro Iniestal Il992l ) for each 
horizontal position and time step as emerging at solar disk center. 

Fig. [1] compares two snapshots of the adiabatic simulations with a vertical 
(left) and horizontal (right) driver. The vertical driver excites in the photosphere 
the fast (acoustic) mode and the surface mode. The fast mode is transformed 
into the slow acoustic mode after the cs = va height and continues up to the 
chromosphere forming shocks with an amplitude of 3-5 km s~^. Note that, in 
the case of the vertical driver, all variations are approximately in phase on both 
sides of the flux tube. The horizontal driver excites in the photosphere the slow 
(magnetic) mode and also the surface mode. The slow mode produces shaking 
of the tube in horizontal direction. After the transformation layer, most of the 
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Figure 2. Time evolution of Stokes V profiles of the Fel 1564.8 nm line in 
the adiabatic simulations with a horizontal (panels a, b, and c) and a vertical 
(panels d, e and f) driver. In each group, the two first panels give spectra 
at two symmetric positions with respect to the tube axis, while the third one 
gives the spatially averaged spectra over the simulation domain. 



energy of the slow magnetic mode goes, again, to the slow acoustic mode. Note 
that, in these simulations with the horizontal driver, the variations are approxi- 
mately in antiphase on both sides of t he flux tube at photospheric heights. For 
more details on these simulations, see iKhomenko fc Colladoi 



^pneric ne 
i fl2007l ). 



3. Surface mode and acoustic mode signatures 

Fig. [2] gives the time evolution of the Stokes V spectra of the Fel 1564.8 nm 
line for different cases (see figure caption). In the case of the simulations with a 
horizontal driver the variations of the synthesized Stokes profiles are due to the 
surface wave. The slow mode is transversal and does not produce any significant 
variations of the vertical velocity. The variations of amplitude and zero-crossing 
shift of the Stokes V profiles are in antiphase. As can be observed from panel (c) 
of Fig. [21 these variations cancel out in the case of observations with insufficient 
spatial resolution. Thus, no time variations of the velocity (nor of any other 
parameter characteristic of the profiles) would be detected in the photosphere at 
the solar disc center if the driver that produces oscillations is purely horizontal. 

In the case of the simulations with the vertical driver, time variations of the 
Stokes V profiles are mainly caused by the fast acoustic mode. The oscillations 
are in phase on both sides of the tube. Thus, spatial averaging does not affect 
the resulting variations to a large extent (Fig. [2]f). An observer would detect low 
spatial resolution velocity oscillations with an amplitude of some 150 m s~^ and 
Stokes V amplitude variations of 10~^, in units of the continuum intensity. Note 
also that, according to Fig. [1] (left panel), the phase of the oscillations changes 
rather slowly with height. Thus, oscillations detected in the photosphere and 
the chromosphere would be coherent (with some delay) if the driver that excites 
oscillations has a vertical component. 



4. Waves at canopies 



Horizontal photospheric shaking of the flux tube due to the slow mode in the 
simulations with a horizontal driver produces strong compression on the inter- 
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Figure 3. Time evolution of the vertical velocity, in a part of the domain 
(300x400 km), from the simulations with a horizontal driver, showing the 
wave traveling at the canopy zone. Numbers at each panel give the elapsed 
time in seconds. The small square refers to a position where velocities are 
similar to those derived from the V-profiles of the Si I line (see Fig. 4). 




Figure 4. Top: diamonds give time variation of the zero-crossing velocity of 
the Si I Stokes V profiles (positive velocities meaning downward motion). The 
solid curve gives the vertical velocity extracted directly from the simulations 
at position marked by black square in Fig. [3l Bottom: time evolution of the 
Stokes V amplitude asymmetry. 



faces between the flux tube and its non-magnetic surroundings. Due to that, 
already at the heights of formation of the Si I line, the vertical velocity varia- 
tions induced by the surface mode show a weak non-linear behaviour, which is 
illustrated in Figs. [3]and[4l 

An inspection of the temporal evolution of the simulations (Figs. [3]) reveals 
an increase of the vertical velocity amplitude just at the border between the flux 
tube and the field free atmosphere and the absence of any significant velocity 
variations outside the flux tube. These abrupt changes produce a weak shock 
wave behaviour of the zero-crossing velocity calculated from the Stokes V pro- 
files of the Si I line (upper panel of Fig. [4]). For comparison, we also show in this 
figure a velocity curve extracted directly from the simulations at the correspond- 
ing height. The similarity between both curves is evident, confirming that the 
variations of the Si I Stokes V zero-crossing position are related to heights about 
300-350 km in the photosphere. The magnetic field discontinuity, together with 
the velocity variations, at the field-free interface produce strong asymmetries of 
the Stokes V profiles. The variation of the Stokes V amplitude asymmetry can 
be as large as ±20% (bottom panel of Fig. [4j). As expected, the asymmetry is 
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Figure 5. Normalized power spectra of the LOS velocity oscillations in the 
Fel 1564.8 nm line (solid line), in the Si I line (dashed line) and in the simu- 
lations at a height of about 1500 km (dotted line), averaged over the central 
part of the tube. Left panel: adiabatic simulations with a vertical driver; 
right panel: simulations with the same driver, but including radiative losses. 



anticorrelated with the zero-crossing velocity, being maximum when the velocity 
reaches maximum upward values. 



5. Power spectra at different heights 

In the above sections, adiabatic simulations have been discussed. The simula- 
tions with a vertical driver including radiative losses are qualitatively similar 
to those without radiative losses. In the photosphere, both the surface and the 
fast modes are generated and fast-to-slow mode transformation is observed at 
the cs = va layer. However, there is an important difference. In the adiabatic 
simulations, the frequency of the oscillations changes with height, due to the 
non-linear generation of the secondary harmonics, from 3 mHz to 6 mHz. The 
oscillations excited by the driver at 3 mHz are evanescent already in the photo- 
sphere and, thus, can not propagate, unlike their higher frequency harmonics. 
This is not the situation if oscillations are allowed to exchange radiative energy 
with their surroundings. Including energy losses is formally equal to reducing 



the e f fective cut-off frequency of the atmosphere by a certain amount (jRobert 



19831 : ICenteno et all [20061 : iKhomenko et al.ll2008l ). Thus, waves at 3 mHz can 



propagate through the whole atmosphere. 

Fig. Ogives the power spectra of oscillations at different heights. The solid 
curve is obtained from the Fel line profiles, i.e. corresponds to the deep photo- 
sphere. There, the dominating frequency is that of the driver, 3 mHz. Higher 
up, in both, adiabatic and non-adiabatic vertically-driven, simulations, power 
at higher frequencies starts to appear. This can be appreciated from the dashed 
curve obtained from the Si I profiles. There is still not much difference between 
the adiabatic and non-adiabatic simulations at photospheric heights. The dot- 
ted curves correspond to the vertical velocity variations at a height of about 
1500 km, in the chromosphere. Now, the adiabatic simulations show a clear 
shift of the power to higher frequencies, while non-adiabatic simulations have 
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the dominant oscillations at the frequency of the driver at 3 mHz. This implies 
that vertical thin magnetic structures, such as flux tubes, can channel 5-minute 
oscillations from the photosphere to the chromosphere if radiative effects are 
taken into account. 

6. Conclusions 

In this paper we have studied the wave behaviour in small-scale flux tubes by 
means of non-linear numerical simulations and Stokes diagnostics. The following 
items summarize our conclusions. 

(i) If the driver that excites oscillations has a vertical component, velocity 
oscillations with an amplitude of some 150 m s~^ and Stokes V amplitude os- 
cillations of 10~^, in units of the continuum intensity, may be detected in low 
spatial resolution observations. The oscillations are coherent through the whole 
atmosphere (with the corresponding time delay due the upward propagation). 

(ii) If the driver that excites oscillations is purely horizontal, no variations 
would be detected in observations with reduced spatial resolution, since the 
vertical velocity variations are produced by the surface mode and are in antiphase 
on both sides of the tube. 

(iii) The LOS velocity shows a weak non-linear behaviour in canopy regions, 
already at heights of 300-350 km in the photosphere. 

(iv) It is important to take into account the radiative losses of oscillations. 
They can produce a decrease of the effective cut-off frequency, making possible 
the propagation of the 5-min oscillations in the chromosphere in vertical small- 
scale magnetic structures. 
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